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ABSTRACT 

We cond uct Smoothed Particle Hyd rodynamics simulations of the 'collect and collapse' 
scenario (Elmcgrc en fc Ladalll977t ) for star formation triggered by an expanding HII 
region. We simulate the evolution of a spherical uniform molecular cloud with an ion- 
ising source at its centre. The gas in the cloud is self-gravitating, although the cloud 
is prevented from globally collapsing. We find that the shell driven by the HII region 
fragments to form numerous self-gravitating objects. We repeat our calculations at 
four numerical resolutions to ensure that they are converged. We compare our results 
to the analytical model of IWhitworth et all (Jl994) and show that our simulations and 
the predictions of Whitworth et al are in goo d agreement in the sens e that the shell 
fragments at the time and radius predicted by IWhitworth et al.l (|1994[ ) to within 20% 
and 25% respectively. Most of the fragments produced in our two highest resolution 
calcul ations are approximately half the mass of those predicted by IWhitworth et al.l 
(|l994f) . but this conclusion is robust against both numerical resolution and the pres- 
ence of random noise (local fluctuations in density of a factor of ~ 2) in the initial 
gas distribution. We conclude that such noise has little impact on the fragmentation 
process. 
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1 INTRODUCTION a uniform gas distribution, the initially-plentiful supply of 

photons drives a spherical ionisation front radially outward 



The question of whether the extreme mechanical and radia- 
tive luminosities of OB stars can trigger or induce the birth 



from the source at a velocity much higher than the speed of 
sound in the cold neutral gas. T his is known as the R-type 



of other stars is a crucial and still open question in the field of , r , , r , , , . . rrr~^ — I J, r.t .k 1^ , , , I 

; 7ttt phase of the front s evolution. (Kahn (1954), Ostcrbrock 

star formation. The winds, expanding HII regions and even- ]., „„ .k x , , c , , . , , , . . ,. r 

° ° (1974)) As the volume ol gas behind the ionisation front 

tual supernova explosions generated by massive stars sweep ,. . , , r 

° . , grows, recombinations consume a progressively farger trac- 

up and compress the surrounding gas m their natal molec- , . , , , , , , . , „ 

f ° tion or the ionising photons and the propagation speed ot 

ular clouds. The Jeans mass, Mj, denning the maximum , , , „„ . , , ., r . , 

° the front drops. When the expansion speed ot the front be- 

mass of gas at a given temperature stable against gravita- , n n n . . , , , „ 

° . , , , comes comparable to the sound speed m the gas, the front 

tional collapse, is related to the gas density n as Mj oc 1 ufp ... r , . „ , . ^ , . , . , , , , , , , 

° f i i transitions from being K-type to D-type m which the bubble 

(for a spherical obicct). Isothermal compression of molecular r . . , , . , , . , 

v . of ionised gas expands at the sound speed m the ionised gas 

gas therefore decreases the Jeans mass, encouraging gravi- /mi -l c i_ j ± ir.4 mi j- ^- ^ 

° ' a ° ° (~ 10 kms for hydrogen at 10 K). The radius at which 

tational fragmentation and, ultimately, star formation. ... . ... . „. .. n . „ . , 

° ' this transition occurs is the atromgren radius, K s , given by 
The collect a n d col lapse model, first proposed by 

lElmegreen fc Ladal (|l977h . is one of the simplest mod- j / L* \ 1/ ' 3 f n. \ ~ 2 / 3 



L ( L* \ ' / n \ * 

R °= Vi^~ a6 H io%i ) liow^J pc - (1) 



els for induced star formation, since it is essentially one- 
dimensional. It therefore provides a useful test-case in which 
theoretical, numerical and observational analyses may be where L t is the source luminosity in ionising photons per 
relatively easily compared. second, n is the initial number density and olb is a recombi- 

When an ionising source is suddenly turned on within nation coefficient, defined such that the recombination rate 

per unit volume to all atomic states except the ground state 
is otBn 2 . Recombinations direct to the ground state, which 
* E-mail: Jim.Dale@astro.le.ac.uk (JED) produce secondary ionising photons, are treated by the 



2 J. E. Dale, I. A. Bonnell, A. P. Whitworth 



'on-the-spot' approximation - secondary ionising photons 
are assumed to be absorbed immediately very close to their 
site of emission. This assumption is achieved by modifying 
the true recombination coefficient a, yielding the modified 
recombination coefficient a_g. The Stromgren radius thus 
simply defines the volume in which all the ionising photons 
emitted by the source are consumed by recombinations 
within the ionised gas, assuming that the gas density is 
unchanged by expansion. 

The expansion of the HII region through the sur- 
rounding neutral gas is highly supersonic and drives a 
shock, sweeping up an increasingly massive and dense 
shell of cool neutral gas - this is the 'collect' phase of 
the 'collect-and-collapse' process in which the HII region 
simply acts like a snowplough. 

If the expansion of the HII region continues for long 
enough, the surface density of the shell increases to the 
point where the shell becomes self-gravitating. The shell is 
then expected to fragment and individual fragments may 
then enter a non-linear collapse phase, possibly forming 
stars. This feature of HII region expansion has been rather 
less well studied, owin g to the non-linear nat ure of the 
evolution at this stage. Iwhitworth et al.l (|l994T ) (hereafter 
usually referred to as W94) studied the evolution of an 
infinite uniform volume of gas in which an ionising radiation 
source is suddenly ignited. They used a two-dimensional 
Jeans analysis to study the evolution of the dense shell 
swept up by the HII region. In the next section, we describe 
W94's model, since this forms the basis of our simulations. 
In Section 3, we describe our numerical methods and 
simulations. We present our results in Section 4 and draw 
our conclusions in Section 5. 



2 A THEORETICAL 

COLLECT-AND-COLLAPSE MODEL 

IWhitworth et all (|l994l ) consider the gravitational stability 
of the shocked shell driven by an HII region expanding 
in a uniform medium of density po- We summarise their 
methodology here - for a more detailed exposition, see their 
Section 5 and Appendix B. 

W94 represent the time evolution of the radius of 
the HII region, and therefore of the shocked shell, as a 
power law given by R(t) = Kt 13 . In a u niform medium , 
K « (JaHii/^f^Rj 1 and p = 4/7 (e.g. ISpitzer! |l97St )1. 
For R(t) >a fewx_R s , nearly all the material with r > R(t) 
is neutral material in the swept-up shell, so the surface 
density E of the shell is ~ poR(t)/3. They consider a 
small circular element of the shell and estimate the inward 
(i.e. azimuthal) acceleration experienced by the fragment, 
which is a function of the fragment's self-gravity and 
internal pressure, and of the velocity divergence due to the 
expansion of the shell. They derive the shortest timescale, 
tfrag on which such an element will begin to condense out 
of the shell and take fragmentation to begin at t — tf rag . 
Using the power-law expression for the expansion of the 
shell, they derive expressions for the radius of the shell 
when fragmentation begins Rf rag , and the mass Mj rag 
of these fragments. For convenience, we reproduce the 
expressions derived by W94 below. In these expressions, 



a. 2 is the sound speed inside the shocked layer in units 
of 0.2 km s _1 , L49 is the source ionising flux in units of 
10 49 photons s _1 , and 113 is the initial gas atomic number 
density in units of I0 3 cm -3 . In these simulations, we 
consider gas which is initially pure molecular hydrogen, 
so the effective atomic number density is simply twice the 
molecular number density. 

r> co 4/11 r 1/11 -6/11 ,„) 

Rfrag ~ 5.8a J Z/49 72 3 pC (2) 

i , a 7/11 r — 1/11 -B/ll-w , s 

tfrag ~ 1.6a J £49 % Myr (3) 

M frag ~ 23a 4 2 0/11 L 4 - 9 1/11 n- 5/11 M Q (4) 

3 NUMERICAL METHODS AND 
SIMULATIONS 

We conduct our simulations using the Smoot hed Parti- 
cle Hy drodynamics code described in detail by iBate et"ail 
\ 19951 ). This code is a hybrid N-body-SPH code in which 
dense, gravitationally-bound and contracting clumps of gas 
particles can be replaced by point masses called sink parti- 
cles. To increase efficiency, particles are integrated on their 
own individual timesteps. We modified the code to allow the 
inclusion of point sources of ionising radiation, which may 
be sink particles or arbitrarily-chosen points in the simu- 
lation volume. Our algorithm defines a 'Stromgren volume' 
around the ionising source and heats the gas within that 
volume to 10 4 K. The particles are first sorted by increasing 
radius from the source. Each particle is then assessed to see 
whether it lies in front of or behind the ionisation front. The 
i particles nearest to the vector between the 'target' particle 
and the sou rce are located using a technique similar to that 
presented in lKessel-Devnet fc Burkertl l|2000l ). The positions 
and densities of these i particles are projected onto this vec- 
tor so that the density pi is known at the positions n along 
the vector. The ionising photon flux per unit solid angle S 
reaching the target particle is then calculated as 

5= 4^ _ Z^ as 2 r (i-i)( r * - r (i-i))- ( 5 ) 

i 

The zeroth particle in the sum is the ionising source, 
whose density is taken to be zero. This approximation is 
equivalent to taking the density variations along this vector 
to represent the radial variations in a spherically-symmetric 
gas distribution, and taking each line-of-sight to be inde- 
pendent, consistent with the on-the-spot approximation. If 
S > 0, so that a non-zero ionising flux is reaching the target 
particle, the particle is ionised and heated. A value of S 
less than zero implies that all the ionising photons emitted 
in the direction of the target particle are absorbed before 
reaching it, so that the particle is beyond the ionisation 
front. 

In the simulations presented here, the Stromgren 
volume is updated every 270 yr (i.e. ~ 10 4 times per simu- 
lation). Increasing or decreasing the ionisation timestep by 
factors of two has no discernible effect on either the results 
or the runtimes of our simulations, so the choice of timestep 
is arbitrary. There is presumably a maximum ionisation 
timestep which, if exceeded, no longer gives consistent 
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results, but we did not investigate this. 

Since simulation of an infinite volume of gas is impos- 
sible, and periodic boundary conditions are inappropriate 
for this problem, we had to construct a finite volume of 
gas in which all the physical phenomena in which we are 
interested could be contained and adequately numerically 
resolved, while still making use of realistic quantities for 
variables of astrophysical importance, e.g. gas density and 
source luminosity. 

We used a finite spherical gas distribution, pressure- 
confined at its outer boundary. The gas is pure molecular 
hydrogen initially. We prevented the cloud from radially 
collapsing by exerting a fictitious force on each particle 
equal in magnitude but opposite in sign to the gravitational 
force on the particle due to all the material interior to the 
particle. We evolved a test cloud for a freefall time without 
an ionising source to confirm that the initial conditions 
were stable and found no detectable particle movement save 
for a thin layer on the edge of the cloud where particles 
moved slightly from their initial positions. The resulting 
change in density over a freefall time was very small and 
not important, since we intended to terminate our simula- 
tions before the shock reached the edge of the cloud. The 
fictitious force was also present during our fragmentation 
calculations. The presence of the force will not affect the 
fragmentation process itself, since the expanding shell is 
confined by thermal pressure from the Hll region and ram 
pressure from the undisturbed gas, not by self-gravity 
for the duration of our simulations. The fragmentation 
process is then essentially two-dimensional, as self-gravity 
is effective only in tangential directions. 

The characteristics of the fragmentation process de- 
pend on three parameters - the isothermal sound speed 
in the neutral gas, the source ionising luminosity and the 
initial gas number density. This results in a large parameter 
space, some of which is not accessible to study numerically 
due either to numerical resolution issues or prohibitively 
long calculation times. We therefore do not intend to 
explore this space thoroughly in this paper. Instead, we 
select a point in the parameter space for which we may 
perform calculations over a range of numerical resolutions 
in a reasonable timeframe an d compare our r e sults with 
the theoretical predictions of IWhitworth et afl (1 19941 ). In 
order to reduce the parameter space from three dimensions 
to two, we choose to simulate molecular clouds consisting 
of pure molecular hydrogen at a temperature of 10K. This 
fixes the isothermal sound speed at 0.2 km s _1 , so that a. 2 = 
unity. We now plot the two dimensional parameter space 
and shade regions which satisfy two numerically-driven 
constraints. The first constraint is that the initial Stromgren 
radius be no less than one tenth of the fragmentation radius, 
ensuring adequate resolution of the HII region. The second 
constraint is that the ratio of the fragment mass to the cloud 
mass must exceed 100 /N par t, where N par t is the number of 
particles in the simulation. This constraint ensures that the 
fragments are adequately resolved by at least one hundre d 
particles, a limi t recom mended by iBate fc Burkertl (|l997t ) . 
iBate fe Burkertl < |l99Tt ) stress that, to resolve gravitational 
fragmentation correctly in SPH, the Jeans mass at which 
fragmentation occurs should be equal to the mass of at least 
twice the mean number of SP H particle neighbours, i.e. 100 
particles. iHubber et al.l l|200fj ) find that even gross violation 
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Figure 1. Plot of the luminosity-number density parameter 
space with the constraint that the initial Stromgren radius be 
at least one-tenth the fragmentation radius. The hatched area is 
forbidden by this constraint and the asterisk marks the parame- 
ters used in our calculations. 

of this condition does not lead to spurious fragmentation, 
although the collapse of marginally-unstable (i.e. slightly 
larger than the Jeans mass) objects is suppressed. Based on 
the predictions of W94, we selected initial c o nditio ns and 
numerical resolutions so that IBate fc Burkertl (|l997T )'s limit 
would be comfortably met. 

We can express R s /Rf rag and Mf rag /M t otai as 

R- s on 4/11 r 10/33 -28/11 /„\ 

8.00.2' L 49 ' n 3 (6) 

-ft/rag 

„ 0.22a 2 2 8/11 L- 4/11 n- 24/11 (7) 

Mtotal 

The fragment mass is a fairly weak constraint in a simu- 
lation where several x 10 particles can be used, so we ignore 
it here. We plot the parameter space with the constraint on 
Rjrag in Figure Q] 

We choose a source luminosity of 10 49 s" 1 , consistent 
with an 07-type star, and an initial neutral gas molecular 
density of 100 cm -3 (corresponding to an atomic number 
density in the HII region of 200 cm -3 ), representative of 
low-density molecular c louds. Making use of the analysis by 
IWhitworth et al.l (| 19941 ) , we then obtain the radius at which 
the shocked shell should fragment as Rf rag — 14 pc, and 
an expected fragment mass of Mf rag ~ 48 Mq. For conve- 
nience, we record the physical parameters of the cloud and 
radiation source and the expected characteristics of the frag- 
mentation process in Table [1] 

In numerical work it is desirable to repeat simulations 
over a range of numerical resolutions (a range of particle 
numbers in the case of SPH) to ensure that convergence 
is achieved and that the physical processes of interest are 
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T 


1n 49 „ — 1 
1U S 


n 


200 cm" 3 




14.6 pc 


M c l ou d 


6.4 x 10 4 M Q 


Rfrag 


13.9 pc 


If rag 


3.3 X 10 6 yr 


Mf ra g 


48 M Q 



Table 1. The physical parameters of the cloud and ionising source 
used in our calculations, and the predictions for the radius and 
time at which the shocked shell should fragment, and the mass of 
the fragments, derived from Equations 2—4. 

being resolved. We therefore conducted a series of simu- 
lations using the physical variables given in Table [TJ with 
particle numbers of 0.37,0.74,1.3 and 2.6 x 10 6 particles. 
Note that the smallest mass that can be resolved in our 
lowest-resolution calculation is ~ 17 Mq, well below the ex- 
pected fragment mass. The 2.6 million-particle simulation 
took about six weeks to run on a 4 x 2.4 GHz SUN v40z 
machine. 



4 RESULTS 

4.1 Numerical convergence 

In its early stages, the evolution of the system under consid- 
eration here is essentially one dimensional, so we first exam- 
ine the growth of the HII region in each of our calculations 
and compare it with the well-known Spi tzer solution for the 
evolution of the HII region radius R(t) |Spitzerl ll978): 

m =R,(i+isgy , ( 8 ) 

where c s is the sound speed in the ionised gas. 

In Figure [2] we plot the evolution of the radius of the 
shocked shell with time in our four calculations. The radius 
of the shell is estimated by taking the average distance 
of the 10 3 densest particles in each simulation from the 
radiation source at r = 0, since these particles should 
reliably reveal the location of the shocked shell outside the 
ionisation front. 

The agreement of the four runs with each other is 
impressive in the early stages of the simulations. The four 
plots are initially very difficult to distinguish, although 
they do diverge slightly as the time at which fragmentation 
occurs (~ 3 x 10 6 yr) is approached. We also note that the 
agreement between our calculations and the analytic solu- 
tion improves somewhat as the particle number increases. 
The evolution of th e shells clearl y agrees well with the 
analytical solution of ISpitzerl l|l978h . 

The fragmentation process itself is rather more difficult 
to analyse. When a quantity of gas in our simulations 
begins non-linear collapse, the gas particles are replaced 
with a po int mass once the i r den sity exceeds an arbitrary 
threshold. IWhitworth et all l|l994h make predictions about 
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Figure 2. Comparison of the evolution of the radius of the 
shocked shell in our four calculations. Lines correspond to particle 
numbers (in units of 10 s particles) as follows: dot-dot-dot-dashed 
line: 0.37, dashed line: 0.74, dot-dashed line: 1.3, solid line: 2.6. 
We also plot the Spitzer solution for this system (dotted line) 

the radius the shocked shell should reach before it frag- 
ments, and the time after ignition of the radiation source at 
which this should occur. We therefore compare the radius 
and time at which the first point mass forms in each of our 
calculations with the predictions of W94 in Figures [3] and U 
The fragmentation radii observed all lie inside the 
radius predicted by W94 and within ~ 25% of the expected 
value. The observed fragmentation radii exhibit a spread 
amongst themselves of ~ 10%. The fragmentation times 
observed are also all less than the theoretical prediction 
by approximately 20%, again with an internal scatter of 
~ 10%, and appear to converge as the particle number 
increases. We conclude that our calculations are adequately 
converged and that they have sufficient resolution to make 
a further study of the fragmentation process feasible. The 
fact that the shells in these simulations appear to fragment 
at earlier times (and therefore smaller radii) than predicted 
by W94 is not surprising, since fragmentation is in reality 
not an instantaneous process. In addition, there is the 
possibility that small quantities of noise in the initial 
particle distribution are likely to cause some regions of the 
shocked shell to fragment at earlier times (and thus smaller 
radii) than is predicted theoretically. We investigate this 
issue in the next section. 



5 FRAGMENTATION 

The particles in all the simulations described above were 
initially distributed on uniform hexagonal-close-packed 
grids. To eliminate the possibility that the initial particle 
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Figure 3. Comparison of the radius at which the first point mass 
is formed in the four simulations (crosses) with the prediction 
made by W94 (dashed line). 



Figure 4. Comparison of the time at which the first point mass is 
formed in the four simulations (crosses) with the prediction made 
by W94 (dashed line). 



distribution has some effect on the fragmentation process, 
we repeated the 1.3 x 10 6 particle simulation with the 
particles distributed randomly. This produces fluctuations 
in gas density of a factor of ~ 2. We compare screenshots 
from these two simulations in Figure [5] We find that 
the fragmentation time and radius for this calculation 
(determined in the same manner as above) are respectively 
2.85 x 10 6 yr and 3.29 x 10 19 cm, very similar values to 
those obtained in the smooth simulation (2.65 x 10 6 yr 
and 3.40 x 10 19 cm respectively). We conclude from these 
results that our results are robust against the presence of 
noise in the gas and that the initial particle distribution 
has no significant effect on the evolution of the system. 

Figure [S] shows that the geometry of the shocked 
shell becomes extremely complex even before the first 
gravitationally-bound fragments form at ~ 2.7 x 10 6 yr. To 
study the process of fragmentation in detail, we devised a 
procedure to identify gravitationally-bound objects within 
our simulations. We first sorted all the SPH particles in order 
of decreasing density. We estimated the size of fragments ap- 
proximately by eye from screenshots of our simulations and 
defined a cutoff radius r cut . We experimented with values of 
r cu t of 0.5, 1.0 and 2.0 x 10 18 cm, finding that the choice had 
little influence on our results. We then descended through 
the sorted list of particles, taking each dense particle to be a 
possible seed for a bound clump. We fetched all the particles 
within r cut of each potential seed and sorted these by radius 
from the seed. We then added these particles to the seed one 
at a time and calculated the thermal, kinetic (in the centre- 
of-mass frame of the group of particles) and gravitational 
energies (E th£rm , E kin and E gra v) of the collection of parti- 
cles. This process was terminated when addition of the next 
particle would have rendered the group of particles unbound 



(i.e. when adding the particle would have have resulted in 
the ratio (Etherm + Ekin)/ \E grav \ rising above a value of 
unity), thereby constructing bound objects around each po- 
tential seed. Particles assigned to clumps were flagged so 
that they could not be used as seeds themselves, and the 
sorted list was descended to an arbitrarily-defined lower 
threshold which was progressively lowered until convergence 
on the number and total mass of bound objects was ob- 
tained for each given snapshot from e ach simulation. In keep - 
ing with the mass resolution limit of iBate fc Burkertl (|l997l ) 
mentioned earlier, clumps consisting of < 100 particles were 
rejected. 

In Figure [6] we plot the increase with time of the total 
quantity of self-gravitating gas in all simulations (including 
the noisy 1.3x 10 6 -particle calculation), and in the 2.6x 10 6 - 
particle simulation for comparison. After a slow start, the 
amount of bound material in all three simulations rises ap- 
proximately linearly with time. The two 1.3 x 10 6 -particle 
calculations and the 2.6 x 10 6 -particle calculation agree very 
well with each other, but the 0.74 and in particular the 
0.37 million-particle calculations both show significantly less 
bound mass towards the ends of the simulations (by factors 
of about two and seven respectively). We observe that, in 
all cases, the amount of bound mass is only a fraction of the 
total mass in the shocked shell. Since fragmentation in both 
these simulations occurs at a radius of ~ 3.5 x 10 19 cm, the 
quantity of mass in the shell when fragmentation begins is 
~ 3 x 10 4 Mq and this mass clearly grows with time. We see 
therefore that only ~ 7 x 10 3 Mq ~ 10% of the material in 
the shell is self-gravitating at the onset of fragmentation (a 
time of ~ 2.8 x 10 6 yr). The density of the shell at the onset 
of fragmentation is pf rag ~ 10 -21 g cm -3 . The freefall time 
within the shell is then (37r/(32Gp /ras )) (1/2) ~ 2.1 x 10 6 
yr. Extrapolation of Figure [6] suggests that the quantity of 
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Figure 5. Comparison of column-density maps (viewed along the z-axis) from three epochs of the 'smooth' (left column) and 'noisy' 
(right column) 1.3 million particle calculations. Each image is 4.5 X 10 19 cm on a side. The noise introduced by randomising the initial 
particle positions is clearly visible in the images in the right column, although the morphology of the shocked shell is very similar in 
both calculations. To a large extent, the shock smoothes out the noise as it expands. 
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Figure 8. The mass functions of bound clumps identified by the method described in the text at six epochs in all simulations. The 
epoch at which each mass function was generated is indicated across the top of the diagram, and the run to which it belongs is given on 
the right. 



bound mass will reach ~ 3 x 10 4 M Q ~ 2 x 10 6 yr after the 
onset of fragmentation, consistent with the freefall time in 
the shell. Due to the fact that the HII region reaches the 
edge of the simulation volume only ~ 0.2 x 10 6 yr after the 
onset of fragmentation, we cannot follow the fragmentation 
process to completion in these calculations. 

In Figure [7] we plot the increase with time of 
the number of self-gravitating clumps (identified by the 
method explaine d abo ve) in all simulations. Given that 
IWhitworth et alj (|1994h estimate the clump mass in these 
simulations to be ~ 50 Mq , we should observe at least ~ 600 
clumps. We see that that the number of clumps rises approx- 
imately linearly with time (except in the lowest resolution 
run, which never forms enough clumps meeting our reso- 
lution limit to draw a conclusion of this nature). The two 
1.3 x 10 6 -particle calculations and the 2.6 x 10 6 -particle cal- 



culation are again in very close agreement with each other 
and approach the numbers of clumps predicted by W94 at 
the ends of the calculations. The 0.74 and particularly the 
0.37 x 10 6 -particle simulations produce fewer bound objects, 
by factors of less than two and about eight respectively, rel- 
ative to the higher-resolution runs. 

In Figure [8] we plot the clump mass functions in all cal- 
culations at six epochs for comparison. The times at which 
the mass functions are taken are given across the top of the 
diagram and the calculation to which they belong is given 
on the right. The vertical dot-dashed lines represent the 
mass resolution limits of each calculation, taken to be the 
mass of 100 particles. In common with Figures [6] and [7] the 
correspondence between the noisy and smooth 1.3 million- 
particle runs and the 2.6 million-particle run is very strong, 
implying both that we have achieved adequate numerical 
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Figure 6. Evolution with time of the quantities of self- 
gravitating material in all simulations (particle number expressed 
in units of 10 6 particles): 0.37 (solid line), 0.74 (long— dashed line), 
1.3 (smooth) (short— dashed line), 1.3 (noisy) (dot— dashed line) 
and 2.6 (dot-dot-dot-dashed line). 



resolution and that noise in the initial conditions has little 
effect on the evolution of the system in question. The last 
two epochs of these runs in particular all show mass func- 
tions decreasing approximately linearly with mass, although 
the gradient is somewhat more negative in the calculation 
with primordial noise, resulting in slightly fewer fragments in 
the range 30 — 50 Mq . By contrast, the mass functions gener- 
ated in the 0.74 million-particle run show more of a peaked 
morphology, with the peak remaining at a roughly constant 
mass of ~ 15 Mq. Finally, the lowest-resolution calculation 
produces so few bound objects meeting the resolution cri- 
terion that the mass-functions generated can hardly be so 
called. This is a result of the relatively poor mass-resolution 
in this calculation. The 1.3 and 2.6 million-particle runs sug- 
gest that the majority of fragments formed in the system 
under investigation (within the timeframe of these simula- 
tions) have masses < 20 Mq and the mass resolution in the 
lowest-resolution run is ~ 17 Mq, so it is not surprising 
that few fragments are observed in the 0.37 million-particle 
calculation. The mean fragment mass observed in the 1.3 
and 2, 6 million-particle runs is ~ 20 Mq, less than half the 
value of ~ 48 Mq predicted by W94, although these runs 
do produce a few objects of this mass. The density in the 
shell at the point of fragmentation, ~ 10 -21 gem -3 implies 
a Jeans mass of ~ 40 Mq , which is consistent with the frag- 
ment masses we observe and with the estimate of W94 to 
within a factor of two. 
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Figure 7. Evolution with time of the number of self-gravitating 
clumps in all simulations (particle number expressed in units 
of 10 B particles): 0.37 (solid line), 0.74 (long-dashed line), 1.3 
(smooth) (short-dashed line), 1.3 (noisy) (dot— dashed line) and 
2.6 (dot— dot— dot— dashed line). 



6 DISCUSSION 

The convergence of our simulations on consistent values 
of tf r a g and Rfrag implies that we can adequately repro- 
duce the gross evolution of this system for our particular 
choice of gas density and source luminosity. In addition, 
the similar evolution of the fragmenting shell in our 
two highest-resolution calculations implies that we have 
sufficient resolution to study the fragmentation process 
itself. Comparison of the high-resolution calculations with 
a third calculation in which we attempted to influence 
fragmentation by introducing random noise into the gas 
distribution shows that our results are robust against 
numerical noise and do not depend on the initial placement 
of the SPH particles. We conclude that we can usefully 
compa re our results to the predictions of IWhitworth et al.l 
1 19941 ). 

Our simulations yield values for tf rag , Rfrag which are 
respectively within ~ 20% and ~ 25% of those derived by 
W94. This implies that the criteria used by W94 to establish 
when and where the shell fragments are approximately cor- 
rect. The fragment masses we observed are a factor of ~ 2.5 
smaller than predicted by W94. This discrepancy is not 
large, particularly when the necessarily approximate nature 
of their analysis is taken into account. However, it was not 
obvious at the outset that fragmentation in our simulations 
should proceed along the lines envisaged by W94. It may 
ha ve been that th e ionisation-front instabilities discussed 
by iGiulianJ (Il979l ) or noise in the initial gas distribution 
had some unforeseen effect on the fragmentation of the 
shell, but neither of these appear to be the case. Our 
simulations and the work of W94 imply that we have a 
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good theoretical understanding of this (admittedly simple) 
problem and that we can move on to the study of more 
complex and realistic systems. 

The most obvious simplification made in this study is 
that the gas is initially uniform and smooth (although we 
showed that small amounts of noise in the gas do not affect 
the results substantially). Several authors have modelled 
the expansion of HII regions in gas distributions with 
power -l aw radial density p rofiles (e.g. iTenorio-Taele et al.l 
jl986t ). iFranco et alj (1 199Cf n relevant to, for example, the 
expulsion of residual gas from proto-globular clusters. This 
work has not been performed in three-dimensions and 
the gravitational fragmentation of the gas has not been 
included. The work presented here implies that both these 
effects can now be included in realistic simulations. 

We have not explored the possibility that the existence 
of well-defined pre-existing clumps in the neutral gas may 
affect the fragmentation process. The passage of the ioni- 
sation front through or past such clumps may cause them 
to collapse, particularly if they are close to gravitational 
instability already. This is a different and more complex 
problem from the one we study here. The many possible 
characteristics of the clumps - e.g. size, mass, density - 
increase the size of the parameter space enormously. We 
therefore leave the study of this problem to a subsequent 
paper. 

We also neglect here the possibility that the stars 
formed by the collapsing fragments will themselves become 
ionising sources. The fragment mass estimated by Whit- 
worth et al of ~ 48 Mq would produce an O-star with an 
ionising luminosity in excess of that of our source and, if 
we could follow the fragmentation process to completion, 
we may obtain hundreds of such objects. 



7 CONCLUSIONS 

We have shown that the analysis presented by W94 cap- 
tures very well the evolution of a shell of self-gravitating 
gas surrounding an expanding HII region. This work and 
that of W94 shows that, from a theoretical standpoint, the 
collect-and-collapse model is a viable model for triggered 
star formation, at least when the gas from which the stars 
are forming is initially roughly uniform. Although the as- 
sumption that_the_g_as_isumfOTm is rather restrictive, the 
results of iDeharveng et al.l (|2005l ) suggest that the collect- 
and-collapse process does occur. They observed a total of 
seventeen HII regions. For two in particular, Sh 104 and 
RCW 79, they find a central ionising star with a roughly 
spherical HII region surrounded by a ring of dense molecu- 
lar gas. In the case of RCW 79, there are clear signs that star 
formation has occurred and continues to occur in the molec- 
ular gas. They caution, however, that there are no signs of 
star formation in progress inside the ring. This would be 
expected if the geometry were truly spherical - Figure [5] 
shows numerous dense fragments whose projected position 
is inside the shocked shell. This may suggest that the ge- 
ometry of these objects is planar and that we observe them 
pe rpendicular t o the plane. This conclusion was also reached 
bv lLang et al.l (2000) in their study of A-Orionis. These au- 
thors further note that models in which the molecular mate- 



rial around A-Ori was initially uniform do not well explain 
the region's present-day morphology. Such a system would 
make an interesting subject for further study. 

Regardless of this consideration, we can conclude that 
feedback from O-type stars can cause the fragmentation of 
even initially smooth and uniform-distributed gas. This is 
induced star formation in the strongest possible sense, since 
a smooth gas distribution will not otherwise fragment and 
form stars. This is an important conclusion, since it implies 
that any cloud of molecular gas can in principle be made to 
form one or more stars under the influence a nearby O-star. 
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